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1 Introduction 



The interacting boson model (IBM) introduced by Arima and Iachello pQ is 
an algebraic model that has its roots in the nuclear shell model. The approx- 
imation of the IBM that only L = and L = 2 nucleon pairs are considered 
and mapped onto s and d bosons gives rise to the group structure U (6). This 
serves as the dynamical algebra of the model, i.e., the Hamiltonian and other 
operators can be expressed in terms of the generators of U(6). Furthermore, 
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the £7(6) group structure leads to the remarkable property that the Hamilto- 
nian is analytically solvable for certain choices of the interaction parameters. 
In spite of its microscopic underpinning in terms of the shell model of the 
atomic nucleus, the IBM can also be linked to a macroscopic interpretation 
of the nucleus by means of the coherent-state formalism [21I3E] . This formal- 
ism allows one to associate an energy surface in the collective quadrupole 
shape parameters j3 and 7 with any IBM Hamiltonian. Hence, the analyt- 
ically solvable limits of the Hamiltonian, the £7(5), 0(6) and SU(3) limits, 
can be linked to a spherical vibrator, a 7-independent rotor, and a prolate 
or oblate deformed rotor, respectively. The evolution of the IBM Hamiltonian 
with varying parameters and the associated energy surface has been stud- 
ied extensively [5 f 6. 7. 8, 9 J. It was shown that the energy surface undergoes a 
first-order quantum phase transition in the passage from £7 (5) to SU (3) and 
a second-order quantum phase transition from £7(5) to 0(6). The concept of 
quantum phase transitions was introduced by Gilmore et al. [TOllTf] in anal- 
ogy with the well-known thermodynamic phase transitions. Quantum phase 
transitions are not driven by the control parameter temperature, however, but 
rather by the parameters of the Hamiltonian describing the quantum system. 
In the IBM in its simplest form the bosons are restricted to the valence 
space but the model can be extended to a configuration-mixed version (IBM- 
CM) [I2|T3] where particle-hole (p-h) excitations across a closed proton or 
neutron shell are incorporated. In certain regions of the nuclear chart these p- 
h excitations descend very low in energy such that they can strongly interact 
with the regular configuration or even become the ground state. Macroscopi- 
cally this is understood as shape coexistence, the coexistence of several minima 
of the energy surface within a very small energy interval. This macroscopic 
information can be extracted from the IBM-CM by calculating the expecta- 
tion value in a coherent state appropriate for configuration mixing [13]. The 
resulting energy surface exhibits a single minimum or several coexisting min- 
ima depending on the IBM-CM parameters. Recently, the energy surface for 
£7(5)-0(6) mixing has been studied [15] and it was shown that the IBM-CM 
in this case gives rise to an extended phase with shape coexistence. 
The aim of this paper is the study of the more general case of £7 [5)-Q(x)Q(x) 
mixing, where Q(x) ( see sec t- [2]) is the quadrupole operator which drives the 
system to deformation and £7(5) is the spherical-vibrator limit of the IBM. 



2 The energy surface for £7(5)— Q(x)Q(x) mixing 

The most compact form of the IBM Hamiltonian, which captures the essential 
physics of the model, is obtained within the consistent-Q formalism [T6] 

H cqi = en d - \k\Q(x) ■ Q(x) , (1) 
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where hd is the <i-boson number operator and Q(x) — (s^d + cfts)^ + x{ ( $d)^ 
the quadrupole operator. For specific choices of the parameters e, k and x-, the 
three symmetry limits of the IBM are obtained: the U(5) limit for k=0, the 
0(6) limit for e = and x — and the SU (3) limit for e = and x — ±v7 /2. 
By calculating the expectation value of the Hamiltonian ([!]) in a normalised 
projective coherent state [2f3f4"] 
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where N denotes the number of valence bosons and (/3, 7) are collective vari- 
ables. If the values of the parameters for the three different IBM limits are 
inserted, it is found that the U (5) limit can be associated with an energy sur- 
face with a spherical minimum, the 0(6) limit with one with a deformed but 
7-independent minimum, and the SU(3) limit with an energy surface which 
has either a prolate (for x — — v^7/2) or an oblate (for x — deformed 
minimum. 

An extended version of the IBM with configuration mixing (IBM-CM) allows 
the simultaneous treatment and mixing of several boson configurations which 
correspond to different particle-hole (p-h) shell-model excitations [T2|T3] . In 
particular, configurations with N, N + 2, iV + 4, . . . bosons are associated with 
Op-Oh, 2p-2h, 4p-4h, . . . excitations, respectively. In case of mixing between 
a 'regular' Op-Oh and a 'intruder' 2p-2h configuration, the Hamiltonian can 
be written as 

H = P^H^Pn + P ] N+2 {H^t 2 + A) P n+2 + V mix , (4) 

where Pjy and Pn+2 are operators projecting onto the iV-boson and (N + 2)- 
boson spaces, respectively. Although the Hamiltonians are formally equivalent, 
the different superscripts in H^ f and H^~ 2 indicate that the parametrisation 
can be configuration dependent. The parameter A is the energy needed to 
excite two particles across a shell gap, corrected for the pairing interaction and 
a monopole effect [17] . Finally, V m i x = w (s^s^ + ss) + 102(0^ • (ft + d ■ d) denotes 
the interaction between the two configurations. This form of the Hamiltonian 
can be easily extended to incorporate configuration mixing with higher-order 
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particle-hole excitations. 

The geometric interpretation of the IBM-CM is obtained by introducing a 
matrix coherent-state method [Hj. In case of mixing between a Op-Oh and a 
2p-2h configuration the energy surface is given by the lowest eigenvalue of the 
matrix 



u,(P) = (iV,A 7 |y mix |iV + 2,/3,7) = y/(N + 2)(N + l) W0 1 + p , (6) 



and E N (ei, \ki\, Xi, P, 7) and E N+2 (e 2 , |k 2 |, X2, P, 7) the expectation values of 
and H^f~ 2 m * ne appropriate projective coherent state. For simplicity's 
sake, wo and W2 are taken equal (wo = W2 = w) such that u becomes (3 inde- 
pendent. 

In the language of catastrophe theory [18J, which can be used to study quali- 
tative changes in the energy surface, (N, ex, \ki\, Xi, £2, l K 2|, X2, A, uj) are called 
control parameters. Because a general study of the energy surface as a function 
of all nine control parameters is totally out of (current) computational reach, 
we focus on mixing between the dynamical symmetries of the IBM which 
can be considered as the benchmarks of the model. In the present paper we 
concentrate on the case of mixing between the vibrational U(5) limit and the 
deformation driving quadrupole term Q(x) ' Qix) which incorporates both the 
0(6) limit and the SU(3) limit for x — and x — iV</2, respectively. In 
a forthcoming paper the case of mixing between deformed configurations will 
be treated. 

The energy surface resulting from the matrix coherent-state method in the 
case of U(5)-Q(x) ■ Q(x) mixing is given by 



E N 



(ei, |/ci|,xi,&7) 



^ +2 (e 2 ,|«: 2 |,X2,/?,7) + A 



(5) 



with 
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e'N - (N + 2)(1 + x 2 ) - j(N + 2)(N + l) x 2 + A' 



e'N - (N + 2) (6 + x 2 ) - 4(JV + 2){N + 1) + 2A' /3 



+ 



+-(JV + 2){N + l)VUx(3 3 cos(3 7 ) - 5{N + 2) + A' 



e'N +(N + 2)(1 + x 2 ) + -(N + 2)(N + l) x 2 - A' 



+ 



e'N + (iV + 2)(6 + x 2 ) + 4(iV + 2)(JV + 1) - 2A' (3 



(N + 2)(iV + l)VUxP cos(37) + 5(iV + 2) - A' 



+uj' 2 (1 + 



2\4 



(7) 



where A' = A/|/c|, e' = g/\k\ and io' = 2u/\k\. We will omit the scaling factor 
\k\/2 from now on as the structural properties only depend on A', e' and to' . 



3 Criticality conditions and Maxwell points 



3. 1 Introduction 



In the following we rely on the ideas of catastrophe theory as discussed exten- 
sively by Gilmore |18j . In the family of energy surfaces £L (/3, 7; e', \i A', 00', N) 
under study, (e', X , A', lu', N) are referred to as the control parameters while 
(J3, 7) are the collective variables. In general, for an arbitrary set of control 
parameters (e' , X o-> A' , u' Q , N ), the energy surface in 7) exhibits isolated 
critical (or equilibrium) points. Isolated critical points are characterised by 
a vanishing gradient of the energy surface (V-EL = 0) and a non-zero de- 
terminant of the stability matrix (det(«S) 7^ 0). The latter matrix is defined 
as 

f d 2 E^ 



S 



dp 2 
d 2 E_ 



d 2 E_ 

df3d<y 
d 2 E_ 



d^df3 d"f 2 

and its eigenvalues determine the stability properties of the energy surface in 
isolated critical points. If all the eigenvalues of S are positive, the isolated 
critical point is a minimum; negative eigenvalues of S indicate a maximum 
whereas positive and negative eigenvalues characterise a saddle point. Since 
the isolated critical points are the extrema or saddle points of the energy sur- 
face E_(/3, 7; e' , Xo, A' , u' , Nq), they determine its global behaviour. For spe- 
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Fig. 1. Illustration of the fact that degenerate critical points organise the qualitative 
behaviour of a family of functions. Panels (a), (b), (c) and (d) show an arbitrary 
function for different sets of the control parameters. In panels (a), (b) and (d) the 
function exhibits only isolated critical points (indicated with a dot) while in panel (c) 
it has a degenerate critical point (indicated with a square). The number of isolated 
critical points or extrema changes when the control parameters pass through (c) 
corresponding to a degenerate critical point. 

cific values of the control parameters the energy surface exhibits points where 
the determinant of the stability matrix vanishes (det(«S)=0). These are the 
degenerate critical points and they are of great importance. Whereas isolated 
critical points organise the qualitative behaviour of a single energy surface, 
degenerate critical points organise the qualitative behaviour of the entire fam- 
ily of energy surfaces E-(/3, 7; e', x, A', u>', N). If the control parameters are 
varied and pass through values where the energy surface exhibits degenerate 
critical point (s), the topology of the surface changes. This can be understood 
intuitively by realising that the topology of an energy surface is determined 
by its isolated critical points. Consequently, if two or more isolated critical 
points merge into a single degenerate one, the topology of the energy surface 
changes. This is illustrated in fig. [TJwhere the evolution of an arbitrary function 
with varying control parameters is shown. In panel (a) the function exhibits 
4 isolated critical points (indicated with a dot). If the control parameters are 
changed, two extrema (2 and 3) move towards each other (panel (b)) until 
they merge into a degenerate critical point (indicated with a square) in panel 
(c). In panel (d) the degenerate critical point has disappeared and only two 
extrema remain. It is clear that (c) with its degenerate critical point separates 
the region where the function exhibits 4 extrema from the region where it has 
only 2. 

Summarising, the degenerate critical points mark out the different regions in 
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the control parameter space where the qualitative properties of the energy 
surface remain unchanged. Hence, they determine specific lines in the phase 
diagram. If two such lines in the phase diagram intersect, the degeneracy of 
the crossing point is higher than the degeneracy of the critical points deter- 
mining the lines in the phase diagram. This crossing is called a triple point. 
In regions of the phase diagram where the energy surface has several minima, 
it is of interest to know which of these is the global minimum and where it 
jumps from one minimum to another (i.e., where two degenerate global min- 
ima occur). The locus of points in control parameter space where this jump 
of the global minimum occurs is called the set of Maxwell points. In the case 
of U(5)-Q(x) ' Q(x) mixing, the Maxwell points are the solutions of 



dp 



= 0, E_(P = (3 )-E4f3 = 0)=Q . (9) 



3.2 Analytical solution in 7) = (0, mr/3) 



In general, the criticality conditions 



dE dE 

— = 0, — =0, det(5) = 0, (10) 



have to be solved numerically. However, an analytical solution can be found. 
Since mixing between spherical U(5) and deformed Q(x) • Q(x) is considered, 
we expect important changes in the energy surface to occur at (3 = which 
is a minimum in the spherical case and a maximum for the deformed case. 
Expanding the energy surface E- of eq. §f§ around (/?,7) = (0,wr/3) (n 
integer), we find 



E- = t 00 + ^t 20 (3 2 + ^t 30 P 3 + ^UoP* + ^hoP 5 + ^ 32 /3 3 7 2 + ^t 60 p e + 
with 



(11) 
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where the notation ( = —A' + 5(N + 2) is used. 



rl 



(12) 



The criticality conditions (10) are automatically fulfilled in the point (/?, 7) = 
(0,mr/3) as the linear terms in /3 and 7 as well as the quadratic term j 2 
are zero. Consequently, the lines determining the phase diagram in control 
parameter space for (/3, 7) = (0, mr/3) are found by requiring a vanishing 
1 term in the Taylor expansion. Hence, if £20 — 0, the stability matrix S 
vanishes identically and we obtain the locus of four-fold degenerate critical 
points (two-fold in (3 and two-fold in 7), 



(iV + 2)(4iV + x 2 )C + VC 2 + ^ 2 



N 



c - \K 2 + v 



rl 



(13) 



In case of C/(5)-0(6) mixing one has x = and the energy surface ^ will 
exhibit a /3 4 behaviour around (/?, 7) = (0, mr/3), if the control variables are 



chosen according to (13). In all other cases the behaviour of the energy surface 
at the critical points ( 13 ) is of dominant /3 3 character. Note that the global 



behaviour of the analytical critical line remains essentially unchanged when the 
intruder configuration changes from a 7-independent rotor to a prolate/oblate 
rotor as x is P ar t °f a positive scaling factor. If ( > or A' < 5(iV + 2), the 
curve converges asymptotically to the lines &' — and e' = (N + 2)(4iV + 
X 2 )/N. If A' = 5(iV + 2), e' c takes on the constant value (N + 2) (4iV + x 2 ) /N. 
If C < or A' > 5(iV + 2), only the asymptote e' = (N + 2)(AN + x 2 )/N 
remains. As long as e' < e' c , a deformed minimum is found whereas the energy 
surface exhibits a spherical minimum for e' > e' c . If the excitation energy of 
the intruder state goes to infinity, A' — > 00, eq. (13) reduces to e' c = 0. 
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In order to find higher-order degenerate critical points, higher-order terms 
are required to vanish. The coefficient £30 vanishes for x = 0, or for uj' = 
and A' > 5(iV + 2), or for A' = 00. For x — the coefficients t 32 , t 50 and 
all higher-order terms with 7 dependence disappear. This can also be seen 
from the energy surface ^ which becomes 7 independent if x = 0. If we 
additionally impose that £40 = 0, we find the triple point 



4iV 2 (iV + 2) + C(iV + l) 4Ny/{N + 2)[AN 2 (N + 2) + ((N + 1 )] 

W 2 ' n + 1 



(14) 



In the (e', u/) plane of the control parameter space, eq. (14) gives the triple 



point where the analytical solution (13) for 7) = (0, nw/3) and the nu- 



merical solution of eq. (10) for 7) 7^ (0,n7r/3) intersect. From the Taylor 
expansion it is seen that the energy surface exhibits a /3 6 behaviour in the 
vicinity of the triple point. 



3.3 Numerical solutions 



In the general case (/3, 7) 7^ (0,n7r/3), the critical points (see eq. (JlOJ) ) and 
the Maxwell points (see eq. ^ ) must be calculated numerically for the energy 
surface E_ in eq. ([7]). To simplify the numerical treatment, 7 can be "frozen" 
to a certain rnr/3. This follows from the fact that the condition dE_/d , y = 
implies necessarily that x — 0, 7 = rnr/3 or (3 = 0. The case (3 = has 
already been treated in the analytical solution (sect. |3.2| and all 7 dependence 
disappears if x — 0. If we choose 7 = rnr/3, d 2 E_/d r y 2 is positive definite as 
long as 7 and the sign of x an d (3 are chosen consistently^ 1 "] Hence, 7 can be 
fixed to mi /3 without loss of generality and the criticality conditions reduce 
to 

dE_ d 2 E_ , . 

w = iw= - (15) 

In the following sections, we choose 7 = 0. 



4 Phase diagrams for U(5)— Q(x) ' Q{x) mixing 



In sect. |3.2 we have shown that x is P ar t of a positive scaling factor in the an- 



alytical solution of the criticality conditions. The criticality conditions cannot 

1 If 7=0, x must be negative and (5 > 0. If 7 = 7r/3, a positive sign for x an< l P 
must be chosen. 
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Fig. 2. Phase diagrams in the case of C/(5)-0(6) mixing for several values of A' and 
for N = 10. The locus of analytical critical points is shown as a full line, that of 
numerical critical points as a dashed line and that of Maxwell critical points as a 
dotted line. The dots represent triple points. The inset figures illustrate the generic 
shape of the potential as a function of (5 in each of the zones of the parameter space. 
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Fig. 3. Zoom on the part near the origin of the phase diagrams shown in fig. [2} The 
locus of analytical critical points is shown as a full line, that of numerical critical 
points as a dashed line and that of Maxwell critical points as a dotted line. 

be solved in general if x is considered as a symbolic parameter but they can 
if a specific value of x is taken. In this section we discuss the two benchmark 
cases, namely C/(5)-0(6) and U(5)-SU(3) mixing. 



4-1 C/(5)-0(6) mixing 



The case of U(5)-0(6) mixing is obtained by choosing x = in eq. (J7|). In 
figs. [2] and [3] different phase diagrams are shown for several values of A'. The 
number of bosons is iV = 10. For A' = we note a deformed region to the 
left of the analytical solution (full line), a spherical region to the right and 
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above the numerical solution (dashed line) and a region of shape coexistence in 
between the two lines. These three regions meet in the triple point, indicated 
by a dot. As A' increases, the slope of the analytical solution switches sign 
and a second triple point in the origin is created. The first triple point moves 
down the analytical curve until it merges with the triple point in the origin 
(e' t ,u' t ) = (0,0) for A' = [AN 2 (N + 2)/(N + 1)] + 5(iV + 2) (see eq. @) and 
eventually disappears. When the slope of the Maxwell curve (dotted line) 
switches sign, a small spherical region for low e' values starts growing around 
the origin. As A' increases further, the two spherical regions approach each 
other. When these two regions coalesce, the region of shape coexistence is split 
in two. The small coexistence region disappears when the two triple points 
merge, while the large coexistence region with a spherical global minimum 
shifts towards higher e' for increasing A'. The evolution of the phase diagrams 
is similar for all N, although the value of A' where changes occur varies slightly 
with N. 



4.2 U(5)-SU(3) mixing 
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Fig. 4. Phase diagrams in the case of U(5)-SU (3) mixing for several values of A' 
and for N = 10. The locus of analytical critical points is shown as a full line, that 
of numerical critical points as a dashed line and that of Maxwell critical points as 
a dotted line. The inset figures illustrate the generic shape of the potential as a 
function of (3 in each of the zones of the parameter space. 

The phase diagrams obtained in the case of U (5)-SU (3) mixing (x = ~V7 /2), 
shown in fig. |4j are very similar to those obtained for U(5)-0(6) mixing. It 
is important to note the following major differences, however. As the non- 
trivial triple point {uj' t 7^ 0) has been proven to occur only for x — 0, the 
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critical (dashed and full) lines separating the different regions never meet in 
a triple point. This has consequences for the occurrence of shape coexistence 
at realistic values of the parameters when the intruder states lie very high in 
excitation energy. In principle, the small region of shape coexistence for low 
e' values will only disappear if A' goes to infinity. Hence, even for very high 
excitation energies of the intruder states, there will always be a region of shape 
coexistence, however small, for realistic values of e' and uj' . This is an essential 
difference with the case of U(5)-0(6) mixing. 



5 Phase transitions 

Similar to the Ehrenfest classification for thermodynamic phase transitions, a 
classification for shape or quantum phase transitions has been proposed [15] . 
The criterion involves the energy of the global minimum E m i n as a function of 
a control parameter. A shape phase transition is called of zeroth order if E m \ n 
changes discontinuously at the critical point. If the first derivative of E m [ n with 
respect to the order parameter or its second derivative is discontinuous at the 
critical point, the shape phase transition is of first or second order, respectively. 
The first-order phase transitions are characterised by mixed phase regimes, 
i.e. regimes where different phases coexist during the transition. Typical for 
second-order phase transitions is the transition from an ordered to a disordered 
phase or vice versa. 

5.1 First- order phase transitions 

It has been shown that the energy surface associated with the IBM (with- 
out configuration mixing) undergoes a first-order shape phase transition in 
the passage from £7(5) to SU(3) [31120] . This first-order shape phase transition 
occurs when passing through the small region of shape coexistence along the 
transition path. However, this type of shape coexistence is different from the 
one occurring in the phase diagrams shown here, as the latter is the result 
of mixing between regular and intruder configurations. Nevertheless, although 
the underlying physics differs, the similarities in the associated topology sug- 
gest that first-order phase transitions should also occur here when passing 
through the zone of shape coexistence. 

In the left panel of fig. [5] is shown the energy of the extrema of the energy 
surface ^ and of the exact ground state of the IBM-CM Hamiltonian (also 
called the local term of the binding energy [21j) as a function of u/ for N = 10 
bosons, A' = 0, e' = 150, k = —0.01 MeV and x — 0? which corresponds to the 
case of £7(5)— 0(6) mixing. The middle panel displays the first derivative with 
respect to uj' of these quantities for the same values of the other control pa- 
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Fig. 5. Left panel: Energy of the extrema of the potential ([7]) and of the exact ground 
state of the IBM-CM Hamiltonian ([I]) as a function of uj' for N = 10 bosons, A' = 0, 
e' = 150, k = —0.01 MeV and x = 0. The full line corresponds to the energy of the 
spherical extremum (minimum or maximum at (3q = 0) and the dashed line to the 
energy of the deformed minimum (/3o 7^ 0). The global minimum is indicated by a 
thick line. The dots show the evolution of the exact ground-state energy and the 
thin vertical dotted line indicates where the Maxwell line is crossed. Middle panel: 
The first derivative with respect to uJ of the energies shown on the left. The same 
conventions apply. Right panel: The (e',u/) phase diagram with the transition path 
followed in left and middle panels. 

rameters. The right panel shows the path followed in the transition. The local 
term of the binding energy gives the absolute energy of the + ground state. To 
obtain the total binding energy in the IBM, one must also include terms in the 
Hamiltonian which depend only on the number of bosons N [TpT] . As these 
terms are constant for a given N, they do not influence the phase diagram and 
a study of the local term of the binding energy is sufficient. In the left and 
middle panels the full (dashed) lines correspond to the spherical (deformed) 
extrema. From the crossing of the full and the dashed line in the left panel 
of fig. |5j it is clear that the global minimum is deformed until the transition 
path crosses the Maxwell line where the spherical extremum becomes lowest. 
Although the coherent-state formalism can be regarded as a variational mean- 
field method, the exact binding energies cannot be compared directly with the 
energy of the global minimum of the energy surface. This is due to the choice 
of the coherent state which does not carry exact angular momentum L = 0. 
Therefore, the coherent state breaks the 0(3) symmetry of the Hamiltonian 
which is respected in its exact diagonalisation [2]. Hence, rather than an exact, 
quantitative comparison of the energies, the purpose of fig. [5] is to reveal qual- 
itative similarities along the transition path. For a quantitative comparison 
with the exact ground-state energies, the energy surface (j7l) must be projected 
onto L = . From fig. [5] it is clear that the exact ground-state energy and the 
energy of the global minimum evolve similarly with changing control param- 
eter u' . The middle panel of fig. [5] illustrates that the derivative of energy of 
the global minimum (thick line) exhibits a discontinuity at the Maxwell line 
where the energy surface thus undergoes a first-order shape phase transition 
Since the derivative of the exact ground-state energy is reasonably close to the 
derivative of the global minimum and changes rapidly in the neighbourhood of 
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Fig. 6. Left panel: Energy of the extrema of the potential (J7| and of the exact 
ground state of the IBM-CM Hamiltonian ([I]) as a function of e' for N = 10 bosons, 
A' = 406, uf = 250, k = -0.01 MeV and x = -y/f/2. The full line corresponds to 
the energy of the spherical extremum (minimum or maximum at (3q = 0) and the 
dashed line to the energy of the deformed minimum (/3o 7^ 0). The global minimum 
is indicated by a thick line. The dots show the evolution of the exact ground-state 
energy and the thin vertical dotted line indicates where the Maxwell line is crossed. 
Middle panel: The first derivative with respect to e' of the energies shown on the 
left. The same conventions apply. Right panel: The (e',u/) phase diagram with the 
transition path followed in left and middle panels. 

the Maxwell point, we may associate this jump in the derivative of the binding 
energy with a first-order shape phase transition. Summarising, in the case of 
£7(5)-0(6) mixing, the nuclear system undergoes a first-order quantum phase 
transition when passing through the line of Maxwell points. 
The same conclusion is reached for U (5)-SU (3) mixing. In the left panel of 
fig. [6] is shown the energy of the extrema of the potential surface ^ and of 
the exact ground state of the IBM-CM Hamiltonian as a function of e' for 
N = 10 bosons, A' = 406, J = 250, « = -0.01 MeV and x = -V7/2, which 
corresponds to the case of U(5)-SU(3) mixing. The middle panel displays the 
first derivative with respect to e' of these quantities for the same values of the 
other control parameters. The right panel shows the path followed in the tran- 
sition. In fig. [6] the same conventions are followed as in fig. [5j in particular, in 
the left and middle panels the full (dashed) lines correspond to the spherical 
(deformed) extrema. Note that the energies and their derivatives are studied 
as a function of e' whereas the varying control parameter in the case of U(5)- 
0(6) mixing was u'. The fact that the energy surface needs to be projected 
on L = for a quantitative comparison with the exact energies to be valid, 
is immediately clear from the left panel of fig. |6j In the deformed region for 
small e', the exact energy is higher than the energy of the deformed minimum. 
This can be understood by realising that the deformation- driving part of the 
coherent state is also the 0(3)-symmetry breaking part. Hence, the need for 
restoring the 0(3) symmetry by means of angular momentum projection is 
largest in the deformed region. Nevertheless, the global behaviour of the exact 
energy and the energy of the global minimum is similar. The discontinuity 
in the derivative of the energy of the latter (see middle panel of fig. [6]) again 
leads to the conclusion that the energy surface undergoes a first-order phase 
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transition when passing through the Maxwell point. Similarly, the slope of the 
derivative of the exact energy changes strongly in the neighbourhood of the 
Maxwell point and can be associated with a first-order shape phase transition 
of the energy surface. Note that the curve for the deformed minimum exhibits 
a gap corresponding to the passage through the narrow region with a single 
spherical minimum in the phase diagram for A' = 406. 

The cases discussed above are specific examples and can be repeated for any 
set of control parameters passing through a region with shape coexistence. 



5.2 Second- order phase transitions 

In the case of only one configuration, it is known that the transition from 
U(5) to SU(3) is characterised by a first-order shape phase transition. The 
transition from £7(5) to 0(6) on the other hand is of second order [3 20] As 
the energy surface associated with the IBM-CM exhibits a similar transition 
in the case of U(5i)— 0(6) mixing, one expects this transition to be of second 
order. Second-order shape phase transition are recognised from the discon- 
tinuity in the second derivative of the energy of the global minimum with 
respect to a control parameter. Unfortunately, an analysis similar to the one 
for first-order shape phase transitions experiences numerical difficulties in the 
neighbourhood of a critical point. Therefore, we use a different method to 
identify this transition. It is possible to recognise second-order phase transi- 
tions by the observation of a power law behaviour of physical quantities when 
passing through a critical point. In general, a power law describes the power 
behaviour of a physical quantity in the neighbourhood of the critical point, 

F~\a c -a\ 5 (16) 

where F is an order parameter (observable) of the system, a is a control pa- 
rameter, a c the value of the control parameter at the critical point and 5 a 
critical exponent. This behaviour is a fingerprint of second-order phase tran- 
sitions and it is a remarkable fact that phase transitions arising in different 
physical systems often possess the same set of critical exponents. This phe- 
nomenon is known as universality. 

For the specific case of phase diagrams in configuration-mixed systems, we 
study the evolution of the deformation /? at the global minimum of the energy 
surface as a function of the control parameters [5j. We do this for U(5)-0(6) 
mixing and calculate the relevant critical exponents. The main reason for the 
choice of this order parameter is that the evolution of (3q can be treated ana- 
lytically. 

The value of (3q at the global minimum results from solving dE_/d(3 = in the 
unknowns (f3,A',e',(j',N). This equation gives rise to the following relation 
between the control and order parameters: 
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4^-(l + (3 2 )e'N{N + 2)[{N + 2)(3 2 - N] 



(1 + (3 2 ) 2 [[e'N + A(N + 2) 2 }(3 2 + N[d - 4(JV + 2)] 
x ([e'N - 4(JV + 2) + C]/? 4 + [e'N + 4N(N + 2) + 2(\(3 2 + () . (17) 

This relation is derived in Appendix |Aj Since the critical points separating 
the spherical and deformed phases are situated on the locus of analytically 



obtained critical points (13), we focus our attention on the point /3 — and 
e' = e' c (see eq. (13)). For £ ^ (and thus e' c 7^ 4 (AT + 2)) u)± is continuous in 
the point j3 — and e' = e' c . Hence, the sign of u;± in (/3 = 0, e' = ej.) remains 
unchanged in a sufficiently small region around (/3 = 0, e' = e(.). In this point 



^4 reduces to 



4C-v/e'(N + 2) 

Because e' = 4 (AT + 2) is the vertical asymptote in the case of U(5)-0(6) 
mixing, e' c < A(N + 2) if C < and e' c > 4(iV + 2) if C > (See sect. |3~2| ). 



Hence, in the neighbourhood of ((3 = 0, e' = e' c ) cj^ is positive. Expanding uj' + 
around (/3 = 0, e' = ej.), we find to lowest order 



UJ, 



dm' + 



P 



+ 



19V 



2 dp 2 



(3 2 + 



f3=0,e'=e' c 

ld 2 uj' 



dm' + 
~d7 



e - e r 



/3=0,e'=e'<. 



j3=0,e'=e' c 



2 de' 2 



<9V 



J \2 



/3=0,e'=e' c 



df3de' 



/3=Q,e'=e' c 



4Cy/e' c {N + 2) _ 2C(iV + 2)[4(iV + 2) + <] _ ^ . 
e' c -4(iV + 2) Je>c(N + 2)[4(N + 2) ' 



4e' c (Af + 2)[4(iV + 2) + e' c ][4iV 2 (iV + 2) + £(N + 1) - e' c AT 21 



NJe' c (N + 2)[A(N + 2)-e' c ] 2 



P 2 



If we invert relation (13) such that uj' c becomes a function of e' c , we obtain 



0/3 



(19) 



4CJe' c (iV + 2) 



(20) 



e / _ 4 (iV + 2) • 

Again, ui' c with the overall plus-sign is positive. If we subsitute the positive uj' c 
in expression (19), it is clear that u' c and the term u)' + ((3 = 0,e' = e' c ) cancel 
and expression (11 91) reduces to 



4e / c (iV + 2)[4(JV + 2) + e;][4A^ 2 (A^ + 2) + ((N + 1) - e' c N 2 ^ 



N^e' c (N + 2)[A(N 



d\2 



] P 2 



2C(AT + 2)[4(AT + 2) + e / c ] 
e' c (N + 2)[A(N + 2)-e' c ] 2 



(e' c -e'). (21) 
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Fig. 7. Power law behaviour at the critical point (upper panels) and at the triple 
point (lower panels). The left panels shows the deformation at minimum, (3q, as a 
function of e' as predicted by the power law (dotted line) and as compared with 
a numerical calculation (full line). The right panels show the energy surface as a 
function of /3, where the dotted line corresponds to e' < e' c or e' t , the full line to 
e' at the critical or triple point, and the dashed line to e' > e' c or e[. The inset 
plots indicate the transition path in the phase diagram. Other control parameters 
are N = 10, A' = 0, and J = 1000 (upper panels) or J = u' t 930.7954 (lower 
panels). 

We conclude that the deformed minimum occurs with deformation 



2e' c [AN 2 (N + 2) + ((N + 1) - e' c iV 2 ] 1 c ; 

for uj' c > uj[ and e < e' c , (22) 



in the neighbourhood of the critical point (e' c ,u' c , A', N). The conditions uj' c > 
uj[ and e' < e' c ensure that the deformation /3q is real. It is clear that in the 
neighbourhood of the critical point (e' c , u' c , A', N), the deformation /3 at the 
deformed minimum exhibits a power law behaviour with a critical exponent 
1/2. 



If £ = then e' c takes the constant value 4(iV + 2), such that u' ± in (17) is 
discontinuous in (/3 = 0, e' = e' c ). From appendix |A| it follows that u>' is un- 
defined in (/3 = 0, e' = e' c ), hence, it is not possible to derive an analytical 
expression for the powerlaw. 

In the upper panels of fig. [7] the deformation at minimum, /3 , is shown for 
N = 10 bosons, u' = 1000 and A' = as e' crosses the critical line of second 
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order. The evolution of (3q m the neighbourhood of the critical point is com- 



pared with a power law with critical exponent 1/2 as obtained in (22) in the 
left panel. The energy surfaces just before and just after crossing the critical 
point as well as at the critical point are shown in the right upper panel The 
inset figure shows the path followed in the phase diagram. The comparison of 



the power law (22) and the exact (3q demonstrates the validity of the former 



in the neighbourhood of the critical point. 



In the triple point itself (see eq. (14)) the coefficient of (3 in the Taylor ex- 



pansion of uj', in (17) vanishes and higher-order terms have to be considered. 
If C 7^ 0, the Taylor expansion around (3 = 0,e' = e[ becomes to lowest order 
in j3 and e' — e[ 

u' = a 00 + a 01 (e'-e[)+a 02 (e'-e[) 2 + a 21 {e'-e' t )(3 2 + a 03 {e'-e' t ) 3 + a m f3\ (23) 



where 



ooo 



aoi 



0()2 



021 



O03 



4N 
N + l 



(N + 2)[AN 2 (N + 2) + C(N + l)] 



2N 3 (N + 2)[8N 2 (N + 2) + ((N + 1)] 
C(N + l) 2 y/(N + 2)[4N 2 (N + 2) + C(N + 1)] ' 
N 5 (N + 2) 2 {3[8N 2 (N + 2) + C(iV + l)] 2 - 64iV 4 (iV + 2)' 

2( 2 (N + 1)3((JV + 2)[AN 2 (N + 2) + ((N + 1)]) 3/2 
A(N + 2) 2 N 2 [8N 2 (N + 2) + C(N + 1)][AN 2 (N + 2) + C(N 



I)] 1 



040 



C 2 (N + l) 2 ((iV + 2)[4iV 2 (iV + 2) + C(iV + I)])' 7 ' 
N 7 (N + 2) 3 

4( 3 (iV + l) 4 ((iV + 2)[4iV 2 (iV + 2) + C(iV + 1)]) 5/2 

x (5[8N 2 (N + 2) + C(iV + l)] 3 - 64iV 4 (iV + 2) 2 [5((iV + 1) + 24iV 2 (iV 
6(iV + 2) 2 [8iV 2 (iV + 2) + ((N + l)][4iV 2 (iV + 2) + ((N + l)] 2 



2)] 



N(({N + 2)[4iV 2 (iV + 2) + C(iV + 1)] 



3/2 



(24) 



For io' = u>[ the left-hand side of eq. (23) cancels with aoo and the resulting 
equation can be solved as a quadratic equation in /3 2 . Keeping only the leading 
term in (e' — e' t ), we derive the following expression for the deformation f3 Q at 
the minimum in the neighbourhood of the triple point: 



«oi\ 
040/ 



1/4 



c'W 4 



N 4 



1/4 



3(iV + 1) 2 [AN 2 (N + 2) + C(iV + 1)] 



(25) 
(26) 
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Thus, at the triple point the critical exponent changes from 1/2 to 1/4. The 
behaviour of (3q at the triple point and its comparison with the power law of 
eq. (26) is shown in the lower panels of fig. [7j Again, the comparison is very 
good. Note that the critical exponents for the order parameter (3q are the same 
as those in Landau theory of tricritical points [22J . 



6 Conclusion 



In the past years, many theoretical and experimental studies have focused on 
the subject of quantum phase transitions in atomic nuclei. The interacting bo- 
son model provides a tractable framework to study quantum phase transitions 
from a theoretical point of view. Because of the algebraic foundations of the 
model, an energy surface is easily constructed and can be studied within the 
framework of catastrophe theory allowing its qualitative study as a function 
of the control parameters. 

In the present work a detailed study of the energy surface associated with 
configuration mixing between a spherical U(5) and a deformed Q(x) 4 Q(x) 
configuration was performed. By expanding the energy surface around (j3, 7) = 
(0,n7r/3) we have derived an analytical solution of the criticality conditions. 
An analytical expression for the triple point was obtained and it was shown 
that it only occurs in the case of Z7(5)-0(6) mixing. For general (3 7^ the 
criticality conditions must be solved for numerically. The same holds for the 
Maxwell points which indicate where the global minimum jumps from one de- 
formation to another. Phase diagrams for the two most symmetrical cases of 
U(5)-0(6) and U(5)-SU(3) mixing have been constructed an discussed. Both 
cases display a large region of shape coexistence for a broad range of excita- 
tion energies of the intruder configuration. For very high excitation energies 
the presence of the triple point in the case of U(5)-0(6) implies the disap- 
pearance of the region of shape coexistence for low e/\n\ whereas this region 
is always present for other deformed intruder configurations (i.e., for x 7^ 0). 
Finally, we have discussed the order of the shape phase transitions. It turns 
out that, generally, the derivative of the energy of the global minimum of 
the energy surface changes discontinuously at the Maxwell line and undergoes 
a first-order shape phase transition. In a numerical calculation the deriva- 
tive of the binding energy follows this behaviour although the discontinuity is 
smoothed out because of finite-size effects. For the transition from a spherical 
to a deformed minimum in the case of U(5)-0(6) mixing, we have shown that 
the deformation f3 of the global minimum exhibits a power-law behaviour in 
the neighbourhood of the critical point and we have given analytical expres- 
sions for the critical exponents. Hence this transition is of second order. 
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A Analytical solution of dE-/d(3 = in the case of U (5)— 0(6) mixing 



In this appendix we derive the expression (17) which results from solving the 



condition dE_/d(3 = 0. Substituting X = 0, A' = -( + 5(N + 2), u/ = 2lu/\k\ 
and 



ai = e'N + A(N + 2) 
a 2 = e'N - AN(N + 2) - 2( 
b 1 = e'N -4{N + 2) +C, 
b 2 = e'N + AN(N + 2) + 2( 



(A.l) 



in the expression for the energy surface ([7]), we find 



AC 



2(1 +/3 2 ) 2 



ai/3 4 + a 2 /3 2 - C 



(A.2) 



Upon a scaling factor |ac|/2 the first derivative dE_/d(3 can be written as 



dE. 



4(3 



d(3 (1 + (3 2 f 
+ C) 2 W 2 (i + /? 2 ) 4 



( ai /3 4 + a 2 (3 2 - C - [(6i/? 4 + hp 2 
4ai/5 3 + 2a 2 f3 



+ 



l+P 2 ) 2 

(M 4 + b 2 (3 2 + C)(46i/? 3 + 2b 2 /3) + Au ,2 (l + (3 2 f(3 
(1 + (3 2 ) 2 [{b^ + b 2 (3 2 + C) 2 + uj' 2 (1 + /3 2 ) 4 ]5 



(A3) 



This expression can be rewritten to 
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dE. 



2(3 



(1 + /?2)3 [(6^4 + ^2 + Q2 + ^2 (1 + £2)4]^ 

[(2 ai - a 2 )/? 2 + (a 2 + 2C)1 \ihf3 4 + b 2 (3 2 + C) 2 + ^ ,2 (1 + (3 2 f 



x 



+(&i/3 4 + & 2 /3 2 + C) (62 - 2h)f3 2 + {-b 2 + 2C) 



(A.4) 



Assuming that (2a\ — a 2 )(3 2 + (a 2 + 2Q 7^ and /3 7^ 0, the condition dE„/d(3 
leads to 



^ 2fl ,2 l4 = (6^ 4 + b 2 (3 2 + C) 2 [(26! - b 2 )$ 2 + (6 2 - 2Q] 2 
W [(2 ai -a 2 )/? 2 + (a 2 + 2C)] 2 

-(M 4 + 6 2 /3 2 + C) 2 , 

This can be rewritten as 



(A.5) 
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(&i/3 4 + 6 2 /3 2 + C) 2 
(l + (3^[(2a 1 -a 2 )(32 + (a 2 + 2()f 
+{b 2 - a 2 - 4C)) ([2(6! + Ql ) - (6 2 + a 2 )]/? 2 + (6 2 + a 2 )) . (A.6) 



([2(6 X - - (6 2 - a 2 )]/? 2 



Inserting the expressions for ai,a 2 ,6i and 6 2 , we find eq. (17): 



4 V /-(1 + /3 2 )eW(AT + 2)[{N + 2)(3 2 - N] 
(1 + (3 2 ) 2 {[e'N + 4{N + 2) 2 }(3 2 + N[e' - 4{N + 2)[ 
x ([e'N - 4(iV + 2) + (\(3 A + [e'N + 4N(N + 2) + 2Q(3 2 + () .(A.7) 



If (3 = 0, the condition dE_/d(3 = is automatically fullfilled. Hence, there is 
always an extremum, either a minimum or a maximum, at /?o = 0. Finally, if 
(2ai - a 2 )f3 2 + (a 2 + 2() = 0, it follows that 



4(N + 2)[(N + 2)(3 2 -N] 
N(l + P 2 ) 



(A.8) 



Inserting this relation in dE_/dj3 = 0, we derive the following expression for 
( from the condition {pi(3 A + b 2 [3 2 + = (see eq. A.4) 



_ A(N + 2)(3 2 [(N + 3)/3 2 -2N} 

^ (i + A 2 ) 2 



(A.9) 



The index i has been added to demonstrate that e i: Pi and Q cannot be varied 
independently. The parameter u' however can take on any value. Hence, for 
a given Q, there is a corresponding for which the deformation fa of the 



21 



extremum of the energy surface remains unchanged when u' is varied. 
Another solution to BE Id (3 = follows from the condition (b 2 — 2bi)/3 2 
(-b 2 + 2C) = (see eq. A.4). The deformation # then equals ±JN/(N + 2) 



+ 



and €{ ( A.8[ ) becomes zero 
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